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ABSTRACT 

This  report  presents  a  surveillance  radar  model  to  support  the  ongoing  Multi- 
Sensor  Data  Fusion  performance  evaluation  study  for  potential  application  to  the 
Canadian  Patrol  Frigate  midlife  upgrade.  This  surveillance  radar  model  takes  into  account 
the  sensor’s  design  parameters  and  external  environmental  effects  such  as  clutter, 
propagation  and  jamming.  To  a  large  extent,  the  model  is  based  on  the  available  literature 
and  on  relevant  experience  and  recent  studies  in  the  field.  In  addition,  the  latest  findings 
regarding  the  dominant  perturbing  effects  affecting  the  radar  detection  performance  are 
included.  The  radar  model  can  be  used  to  generate  contacts  and  false  alarms  in  scenarios 
for  multi-sensor  data  fusion  studies  while  the  scenario  is  running. 


RESUME 

Ce  document  presente  un  modele  de  radar  de  surveillance  pour  appuyer  une  etude 
en  cours  portant  sur  revaluation  de  la  performance  de  la  fusion  de  donnees  multi-capteurs 
pour  une  application  eventuelle  pour  la  modernisation  de  la  Fregate  de  patrouille 
canadienne.  Le  modele  tient  compte  des  parametres  du  capteur  proprement  dit  et  des 
effets  de  1’environment  tels  le  fouillis  de  mer,  la  propagation  et  le  brouillage.  Le  modele 
s’appuie  en  large  partie  sur  la  documentation  disponible,  mais  incorpore  les  r6centes 
connaissances  relatives  aux  effets  les  plus  perturbants  affectant  la  detection  radar.  Le 
modele  est  approprie  pour  des  dtudes  en  fusion  de  donnees  qui  necessistent  la  generation 
de  vrais  contacts  et  de  fausses  alertes  durant  Texecution  d’un  scenario  donne. 
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EXECUTIVE  SUMMARY 

A  major  ongoing  activity  undertaken  by  the  Data  Fusion  Group  in  the  Command 
and  Control  Division  at  Defence  Research  Establishment  Valcartier  (DREV)  is  to  explore 
sensor  management,  integration,  and  data  fusion  concepts  that  could  be  applied  to  the 
current  Canadian  Patrol  Frigate  (CPF)  Above  Water  Warfare  (AWW)  sensor  suite,  as  well 
as  its  possible  future  upgrades,  in  order  to  improve  its  performance  against  current  and 
future  threats.  The  approach  retained  for  this  activity  is  to  employ  sufficiently 
representative  sensor  and  phenomenological  simulators. 

This  report  presents  a  surveillance  radar  model  to  support  an  ongoing  Multi- 
Sensor  Data  Fusion  (MSDF)  performance  evaluation  study  for  potential  applications  to 
the  Canadian  Patrol  Frigate  midlife  upgrade.  This  surveillance  radar  model  takes  into 
account  the  sensor’s  design  parameters  and  external  environmental  effects  such  as  clutter, 
propagation  and  jamming.  The  latest  findings  regarding  the  dominant  perturbing  effects 
affecting  the  radar  detection  performance  are  included.  In  addition,  the  radar  model  can  be 
used  to  generate  contacts  and  false  alarms  in  scenarios  for  multi-sensor  data  fusion  studies 
while  the  scenario  is  running. 

The  radar  simulation  provides  a  “realistic”  representation  of  a  radar  behavior  in 
typical  operational  environments.  A  panoramic  pulse-MTI/Doppler  search  radar,  typical  of 
the  Canadian  Navy  surveillance  radars,  using  a  rotating  antenna  that  scans  360  degrees  in 
azimuth  is  considered.  The  scenarios  of  interest  consist  of  point-target  (missiles)  attacks  in 
an  open-sea  environment.  The  proposed  model  is  more  versatile  and  more  generally 
applicable  than  those  published  so  far;  it  considers  a  wider  range  of  parameters  which 
affect  a  search  or  surveillance  radar’s  performance.  It  is  therefore  likely  to  yield  more 
realistic  results. 

The  results  of  this  research  are  of  prime  importance  when  considering  mid-life 
update  of  the  Canadian  Patrol  Frigate. 
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1.0  INTRODUCTION 

A  major  ongoing  activity  undertaken  by  the  Data  Fusion  Group  in  the  Command  and 
Control  Division  at  Defence  Research  Establishment  Valcartier  (DREV)  is  to  explore 
sensor  management,  integration,  and  data  fusion  concepts  that  could  apply  to  the  current 
Canadian  Patrol  Frigate  (CPF)  Above  Water  Warfare  (AWW)  sensor  suite,  as  well  as  its 
possible  future  upgrades,  in  order  to  improve  its  performance  against  current  and  future 
threats.  The  approach  retained  for  the  activity  is  to  employ  sufficiently  representative 
sensor  and  phenomenological  simulations.  To  achieve  this  capability,  a  highly  modular, 
structured,  and  flexible  simulation  environment  has  been  developed. 

DREV’s  Concept  Analysis  and  Simulation  Environment  for  Automatic  Target 
Tracking  and  Identification  (CASE_ATTI)  (Ref.  1)  provides  the  algorithm-level  test  and 
replacement  capability  required  to  study  and  compare  the  applicability  and  performance  of 
advanced,  state-of-the-art  Multi-sensor  Data  Fusion  (MSDF)  techniques.  The 
CASE_ATTI  system  allows  the  user  to  create  and  edit  multi-sensor  multi-target  test 
scenarios.  A  typical  scenario  consists  of  platforms  (e.g.,  ships)  with  sensors,  and  targets. 
The  platforms  can  be  stationary  or  moving  along  predefined  paths.  One  or  several  sensors 
(potentially  dissimilar  such  as  radar,  infrared,  ESM,  etc.)  can  be  assigned  to  each  platform. 
Targets  are  created  with  defined  trajectories  and  attributes.  The  scenario  also  allows  for 
the  selection  of  environmental  conditions  that  may  affect  the  various  sensor  detection 
performance  and  measurements  accuracy. 

The  CASE_ATTI  tool  comprises  three  major  modules,  as  shown  in  Fig.l.  The 
sensor  module  is  responsible  for  providing  realistic  measurement  data  to  the  tracking 
module.  This  module  also  offers  the  ability  to  integrate  real  sensor  data  into  the 
simulation.  The  tracking  module  supports  a  wide  variety  of  fusion  architectures,  varying 
from  a  simple  single  sensor  tracker  to  an  arbitrarily  complex  hierarchical  multiple-sensor 
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topology.  The  graphics  module  provides  the  tools  to  examine  the  tracking  performance 
and  the  results  of  the  sensor  fusion.  The  object-oriented  design  allows  the  users  to  easily 
develop  and  incorporate  their  own  tracking  algorithms,  sensor  models,  and  analysis  tools. 

This  report  presents  a  surveillance  radar  model  for  the  CASE_ATTI  sensor  module. 
A  thorough  demonstration  of  its  validity  under  various  threat  and  environmental 
conditions  will  be  presented  shortly  in  a  follow-up  report,  emphasizing  the  case  of  current 
Canadian  Navy  radars.  This  surveillance  radar  model  takes  into  account  the  sensor’s 
design  parameters  and  external  environmental  effects  such  as  clutter,  propagation 
(multipath  and  refraction)  and  jamming.  To  a  large  extent  the  model  is  based  on  the 
available  literature  and  recent  studies  in  the  field.  In  addition,  the  latest  findings  regarding 
the  dominant  perturbing  effects  affecting  the  radar  detection  performance  are  included  in 
the  simulation. 

A  requirement  for  data  fusion  studies  is  that  the  radar  model  can  be  used  to 
generate  contacts  and  false  alarms  in  scenarios  for  multi-sensor  data  fusion  studies  while 
the  scenario  is  running.  Recent  surveillance  radar  simulations,  such  as  CARPET  (Ref.  2) 
and  previous  ones  (Refs.  3-4),  are  not  fully  suitable  for  sensor  fusion  studies  mainly 
because  of  this  requirement.  CARPET  only  gives  plots  of  radar  detection  probability 
versus  range  as  well  as  plots  of  intermediate  results  (or  dependent  parameters),  such  as 
MTI  (Moving  Target  Indicator),  improvement  factor  and  propagation  loss.  The  results 
provided  by  CARPET  are  not  available  as  numerical  data  and  this  prevents  their  use  in  the 
CASE_ATTI  tool. 

The  model  proposed  here  is  made  of  two  levels  of  processing:  off-line  and  on-line. 
At  the  first  level,  various  software  modules  are  provided  to  precompute  tables  that 
represent  the  dominant  perturbing  effects  on  the  radar  detection.  At  the  second  level, 
which  resides  in  CASE_ATTI,  the  probability  of  detection  (Pd  )  is  evaluated  using  the 
tables  each  time  a  target  meets  the  radar  beam.  Measurements  are  then  generated 
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according  to  the  Pd  and  the  signal-to-interference  ratio.  In  addition,  the  simulation 
generates  false  alarms  as  caused  by  clutter  and  system  noise. 

The  radar  model  provides  a  realistic  representation  of  a  radar  behavior  in  typical 
operational  environments.  A  panoramic  pulse-MTI/Doppler  search  radar,  typical  of 
Canadian  Navy  surveillance  radars,  using  a  rotating  antenna  that  scans  360  degrees  in 
azimuth  is  considered.  The  scenarios  of  interest  consist  of  point-targets  (missiles)  attacks 
in  an  open-sea  environment.  The  proposed  model  is  more  flexible  and  more  generally 
applicable  than  those  (Refs.  2-4)  published  so  far;  it  considers  a  wider  range  of  parameters 
which  affect  search  or  surveillance  radar  performances.  It  is  therefore  likely  to  yield  more 
realistic  results.  In  particular,  this  model  brings  the  following  features: 

•  extension  of  sea-clutter  models  to  account  for  ducting; 

•  better  identification  of  the  various  losses  of  typical  surveillance  radars; 

•  full  Parabolic  Equation  Method  (PEM)  computation  of  the  propagation  factor; 

•  a  model  to  generate  false  alarms  induced  by  sea  clutter; 

•  a  flexible  and  generic  model  to  represent  the  effects  of  signal  processing; 

•  a  simulation  architecture  allowing  refinement  in  modeling  as  well  as  real-time  use 
for  data  fusion  studies. 

This  document  is  organized  as  follows.  Chapter  2  gives  a  general  description  of  the 
radar  model.  Chapter  3  presents  the  calculations  of  the  target  signal  power  and  the 
receiver  noise.  In  Chapter  4,  the  interfering  signals  from  clutter  and  jamming  are 
discussed.  Chapter  5  addresses  the  radar  signal  processing  issues.  Chapter  6  and  7  show 
how  the  target  contacts  and  false  alarms  are  generated.  Finally,  Chapter  8  ends  up  with 
the  conclusion. 

This  work  was  carried  out  at  DREV  between  1992  and  1994  under  PSC  12C,  Ship 
Combat  System  Integration. 


Simulation  Manager 


Advanced 
User  Interface 
System 


FIGURE  1  -  CASE_ATTI  global 
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2.0  GENERAL  DESCRIPTION  OF  THE  RADAR  MODEL 

The  search  process  consists  in  scanning  a  volume  of  space  and  reporting  detections 
within  the  volume.  Detections  are  declared  when  the  signal  received  and  processed 
exceeds  a  threshold.  These  detections  can  be  due  to  the  presence  of  real  targets  (true 
detections)  in  the  volume  or  due  to  noise  or  returns  from  objects  that  are  not  targets  (false 
alarms). 

A  panoramic  pulse-MTI/Doppler  search  radar,  typical  of  the  Canadian  Navy 
surveillance  radars,  using  a  rotating  antenna  that  scans  360  degrees  in  azimuth  is 
considered.  The  operation  of  a  typical  pulse  radar  receiver  can  be  described  with  the  aid  of 
the  block  diagram  shown  in  Fig.  2. 


Detection 


FIGURE  2  -  Block  diagram  of  a  pulse  radar 


A  single  antenna  (and  its  lines/waveguides)  is  generally  used  for  both  transmitting 
and  receiving  via  a  duplexer.  The  receiver  is  usually  made  of  a  low-noise  RF  amplifier 
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an  Intermediate  Frequency  (IF).  The  IF  amplifier  is  designed  as  a  matched  filter  to 
maximize  the  signal-to-noise  ratio.  The  IF  signal  is  further  down  converted  to  video  where 
signal  processing  can  be  more  conveniently  applied. 


The  surveillance  radar  model  takes  into  account  not  only  the  sensor's  design 
parameters  but  also  the  determinant  factors  affecting  detection  such  as  clutter,  multipath, 
ducting,  and  the  signal  processing  gain  of  the  radar-environment-target  chain  by 
considering  a  Signal-to-Interference  ratio  (S/I): 


S  S 
I  ~  N0  +  C+J 


[1-] 


where  N0,  C,  and  J  are  the  energy  from  system  noise,  clutter  and  jamming,  respectively. 
In  the  simulation,  the  S/I  is  approximated  by  finding  the  dominant  interference:  N0,  C  or  J. 
In  our  model,  when  these  three  sources  of  interference  are  competing,  clutter  is  taken  as 
being  the  dominant  one. 

The  S/I  is  used  to  evaluate  the  probability  of  detection,  Pd,  each  time  a  target  is 
geometrically  hit  by  the  radar  beam.  Measurements  of  range  and  azimuth  of  each  target 
are  then  generated  according  to  Pd  and  S/I.  In  addition,  the  simulation  generates 
measurements  (range  and  azimuth)  corresponding  to  false  alarms  caused  by  clutter  and 
system  noise. 

The  proposed  model  (Fig.  3)  has  two  levels  of  processing.  At  the  first  level, 
called  off-line  processing  level,  software  modules  have  been  developed  to  precompute 
tables  giving  the  critical  parameters  required  to  calculate  the  S/I.  These  parameters  require 
extensive  computations;  furthermore,  these  external  modules  can  be  conveniently 
upgraded  as  improved  models  become  available.  At  the  second  level,  called  on-line 
processing  level,  which  resides  within  the  CASE_ATTI  sensor  module,  the  S/I  is 
computed  and  the  probability  of  detection  (Pd  )  is  evaluated  each  time  a  target  meets  the 
radar  beam.  The  two-level  processing  architecture  provides  the  flexibility  of  using  various 
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non-standardized  pieces  of  software  to  precompute  sensitive  parameters  without 
increasing  the  complexity  and  the  execution  speed  of  CASE_ATTI. 


FIGURE  3  -  Block  diagram  of  the  surveillance  radar  simulation 

2.1  Off-Line  Processing  Level 

For  various  sets  of  weather  conditions  and  radar  characteristics,  tables  are 
computed  outside  CASE_ATTI  to  provide  parameters  or  coefficients  for  the  computation 
of: 

-  sea-clutter  power  and  skewness; 

-  weather-clutter  power; 

-  propagation  factor; 
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-  filters  improvement  factor; 

-  clutter-induced  probability  of  false  alarms. 

Coefficients  for  sea-  and  weather-clutter  reflectivities  are  calculated  versus  range 
and  azimuth.  For  sea  clutter,  classical  models  have  been  extended  to  incorporate  ducting 
effects.  The  extensions  are  based  on  data  recently  published.  In  addition,  the  skewness  of 
the  statistical  distribution  of  the  sea  clutter  is  estimated  and  described  by  a  factor  and  the 
probability  of  false  alarms  induced  by  sea  clutter  is  also  estimated.  The  propagation 
factor  represents  the  interference  pattern  as  due  to  multipath,  taking  into  account 
refraction  effects.  It  is  computed  using  the  Parabolic  Equation  Method  (PEM)  (Ref.  5). 
The  filter  improvement  factor  is  used  to  compute  the  efficiency  of  clutter  rejection  of 
Moving  Target  Indicator  (MTI)  and  Doppler  filters.  Tables  of  the  filter  improvement 
factors  versus  the  standard  deviation  of  the  Gaussian  clutter  spectrum  and  the  target 
Doppler  shift  are  produced.  The  standard  deviation  and  the  target  Doppler  shift  are  both 
normalized  to  the  Pulse  Repetition  Frequency  (fr ).  Various  kind  of  Doppler  filtering 
processes  are  supported:  iV-pulse  canceler,  TV-pulse  Doppler-filter  bank,  cascade  of  MTI  to 
a  Doppler  filter-bank.  Filter  sidelobe  reduction  can  also  be  provided  by  applying  Taylor, 
Chebyshev,  Hamming  or  Hann  windowing. 

2.2  On-Line  Processing  Level 

In  the  CASE_ATTI  sensor  module,  the  precomputed  parameters  are  used  to 
evaluate  Pd  each  time  a  target  is  geometrically  hit  by  the  radar  beam  and  to  generate  false 
detections.  Of  course,  to  run  a  scenario  in  CASE_ATTI  the  required  precomputed  tables 
for  the  selected  sensors  and  environmental  conditions  must  be  available. 
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2.2.1  Target  Detection 

Each  time  the  beam  hits  a  target,  the  target’s  probability  of  detection  Pd  is 
computed  taking  into  account  the  radar  parameters  and  the  environmental  effects. 
Calculations  are  made  for  the  resolution  cell  that  contains  the  target.  If  more  than  one 
target  happen  to  be  in  the  same  "resolution"  cell,  the  program  checks  if  the  targets  can  be 
differentiated  by  the  Doppler  filters.  The  signal  power  received  from  all  non-differentiated 
targets  is  added.  The  result  constitutes  "the  desired  target  signal". 

Conceptually,  the  S/I  at  the  output  of  the  video  stage  has  to  be  computed,  and  at 
the  detector  module  (thresholding),  the  corresponding  Pd  has  to  be  evaluated.  When  a 
detection  is  declared,  the  radar  model  outputs  the  following  parameters: 

-  target  range  (km); 

-  target  azimuth  (deg-true  north); 

-  signal-to-interference  ratio,  S/I  (dB) 

2.2.2  False  Alarm  Generation 

In  any  radar  system,  a  device  called  the  CFAR  (Constant  False  Alarm  Rate)  is  used 
to  maintain  the  prescribed  false  alarm  rate  constant.  This  device  comprises  a  processor 
which  samples  and  processes  the  interfering  signal  in  the  radar  cells  to  estimate  the  noise 
level  and  dynamically  adjust  the  detection  threshold  to  maintain  a  fixed  Pfa  (probability  of 
false  alarm)  based  on  estimates  of  interference  statistics. 

Noise-induced  false  alarms  are  generated  randomly  at  a  rate  corresponding  to  the 
false  alarm  rate  given  by  the  user.  The  relationship  between  Pfa  and  the  false  alarm  rate 
( rfa )  depends  on  the  radar  characteristics.  In  addition,  further  false  detections  can  be 
generated  due  to  the  spiky  nature  of  sea  clutter.  Most  radar  simulations  have  exclusively 
regarded  the  evaluation  of  target  probability  of  detection  in  noise/clutter,  and  the  current 
knowledge  of  clutter  is  not  sufficiently  advanced  to  provide  an  accurate  representation  of 
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clutter-produced  false  alarm  outcomes.  We  choose  to  use  heuristic  methods  to  achieve  a 
likely  spreading  of  false  alarms  as  a  function  of  radar  characteristics  and  the  environment. 
A  probability  of  false  alarms,  Pfac  is  estimated  as  a  function  of  range  and  azimuth.  Details 
on  the  technique  are  given  in  Chapter  7. 
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3.0  TARGET  SIGNAL  POWER  AND  RECEIVER  NOISE 


This  section  aims  at  computing  the  signal-to-noise  ratio  at  the  output  of  the  IF 

stage. 

3.1  Target  Signal  Power 

The  power  received  at  the  input  of  the  RF  stage  (Fig.  2)  from  a  single-target  echo 
is  given  by 


[2-] 


with 


K  ~  P&K 


(4*) %„ A 


2 

\ atm 


[3.] 


where 

P transmitted  pulse  power  (w), 

G:  maximum  antenna  gain, 

F :  propagation  factor, 
o:  target  cross  section  (m2), 

Lt/r :  transmit/receive  losses, 

Latm :  atmospheric  loss  (oxygen,  water  vapor  absorption  and  rain  attenuation), 

R:  target  range  (m), 

A,:  wavelength  (m). 

3.1.1  Transmit/Receive  Losses 

The  total  transmit/receive  system  loss,  LUr,  is  the  product  of  numerous  individual 
loss  factors  from  the  output  of  the  transmitter  to  the  antenna.  Among  them  are  the  loss  in 
transmission  lines,  the  radome  loss,  the  loss  in  connectors  and  the  attenuation  when  the 
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signal  passes  the  duplexer  (insertion  loss).  It  is  the  responsibility  of  the  user  to  estimate  the 
various  contributions  and  to  input  the  resulting  loss  in  the  radar  specifications. 

3.1.2  Propagation  Factor 

The  propagation  factor  F  accounts  for  the  radar  wave  propagation  effects  resulting 
from  non-free-space  conditions.  It  includes  the  combined  effects  of  refraction  (bending  of 
rays),  multipath,  as  due  to  earth  reflections,  and  antenna  gain  directivity  g( (j>)  (discussed  in 
Section  3.1.3).  The  basic  definition  of  the  propagation  factor  is: 


where  E  is  the  actual  intensity  of  the  electric  field  impinging  on  a  target  at  distance  R  from 
the  source,  while  Eo  is  the  field  that  would  exist  at  the  target  range  in  free  space  with  the 
target  in  the  direction  of  the  maximum  antenna  gain. 

To  perform  field  strength  calculations,  the  use  of  the  parabolic  approximation  of 
the  elliptic  scalar  wave  Helmoltz  equation,  commonly  known  as  Parabolic  Equation 
Method  (PEM),  has  received  wide  acceptance  in  recent  years  (Refs.  5-6).  The  PEM 
approach  applies  a  full-wave,  forward-scatter  model  capable  of  predicting  propagation  for 
arbitrary  atmospheric  refractivities.  In  the  PEM,  the  field  is  expressed  in  terms  of  a 
complex  electrical  field  attenuation  function  in  space,  u(R„  h )  (V/m1/2): 

^2  -v 

^+2yK^  +  K2(m2(R,,/0-l)K  =  0  [5.] 

ah  oR, 

where  k  is  the  free-space  wave  number,  Rt  is  the  terrestrial  range  and  h  is  the  height 
above  the  surface.  The  modified  refractive  index,  m,  accounts  for  the  earth's  curvature  (m 
=  n  +  h/r )  where  n  is  the  refractive  index  and  r  the  earth  radius.  The  propagation  factor  is 
related  to  u  by: 


UNCLASSIFIED 

12 


P154344.PDF  [Page:  29  of  86] 


UNCLASSIFIED 

13 


F  =  \u(R"h)\Jli; 


[6.] 


A  recently  developed  implementation  of  the  split-step  Fourier  technique  is  used  to 
solve  [5]  (Ref.  6).  The  calculation  requires  as  inputs:  radar  frequency  and  polarization, 
antenna  height  (above  the  sea  surface),  antenna  directivity  pattern,  sea-wave  height  and 
refractive  index  of  air.  Arbitrary  vertical  profiles  of  m  can  be  user-defined.  Otherwise,  the 
following  common  default  profile  is  used: 


m(h)  —  1  + 


xHT* 


[7-] 


where  the  evaporation  duct  height  is  hduct{ m).  Mo  is  the  surface  modified  refractivity  (made 
equal  to  320  M  units)  andzo  is  the  roughness  length  equal  to  1.5  x  lO^m.  The  duct  height 
is  calculated  from  the  air  temperature  and  humidity,  the  wind  speed  and  the  water 
temperature  using  the  marine  boundaiy  layer  model  Wamsley-KEL-DREV  (WKD) 
(Ref.  7). 

Figure  4  shows  the  effects  of  ducting  on  the  propagation  factor  against  range. 
Cases  of  duct  heights  of  0,  5,  10  and  15  m  are  shown.  The  target  height  and  the  antenna 
height  are  25  m  and  20  m,  respectively.  We  observe  that  nulls  move  right  as  a  function  of 
duct  height  and  the  null  displacement  becomes  significant  at  ranges  as  close  as  5  km. 

As  mentioned  in  Chapter  2,  the  propagation  factor  values  are  precomputed  in  the 
space  domain  (Ru  h)  for  specific  radar  and  environmental  parameter  inputs.  In  a 
simulation  run,  when  a  table  is  missing  or  a  requested  coordinate  is  located  outside  the  (Rt, 
h)  domain  of  the  supplied  table,  a  standard  multipath  factor  calculation  is  performed,  as 
detailed  in  Ref.  8.  It  is  the  responsibility  of  the  simulation-user  to  provide  the  necessary 
table. 
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3.1.3  Antenna  Gain  Directivity  Factor 

Radar  antennas  usually  have  highly  directional  properties,  with  the  radiation 
concentrated  in  a  beam  in  a  particular  direction.  The  effective  antenna  gain  is  a  function  of 
target  azimuth  0  and  elevation  (j>  and  is  described  by  a  directivity  factor,  g,  of  the  antenna 
pattern.  Directivity  factors  are  required  in  the  calculation  of  the  propagation  factor  and  the 
clutter  and  jamming  signals.  Using  cp  as  a  parameter  that  can  represent  either  the  azimuth 
0  or  the  elevation  <j),  four  typical  directivity  factors  have  been  programmed:  Isotropic, 
Gaussian,  Sinc((p),  and  Cosec. 

The  first  pattern  is  omnidirectional,  the  directivity  factor  is  given  by 

g(<p)  =  l  [8-] 

The  second  pattern  is  based  on  the  assumption  of  a  Gaussian  beam  shape;  it  is 
normally  used  to  describe  a  pencil  beam 

g(cp)  =  exp^-4^2-111  a/2"J  [9.] 


where  cpb  is  the  half-power  beamwidth. 

The  third  pattern  is  a  sine  function;  it  is  used  to  describe  a  pencil  beam  when 
sidelobes  are  needed: 


with 


g(9)  = 


1  +  cos  cp  sin  ( k  sin  cp) 
A:  sin  cp 


[10.] 


A:  s  1.3916/  sin(cpft  /  2) 


[11.]. 


This  pattern  can  be  modified  to  simulate  actual  patterns  having  sidelobe  levels 
considerably  lower  than  13.26  dB  below  the  main-beam  power  level  by  applying  a 
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sidelobe  reduction  factor.  The  antenna  gain  directivity  factor  is  defined  inside  the  main 
beam  by  [10]  and  at  angles  greater  than  the  first  null  which  occurs  at 

tpnu„  =  sin -’(2. 2575  sin  (cp,/ 2))  [12.] 

the  directivity  factor  is  given  by 

g(cp)  =  f1  +  C0S(P^  sin  sin  <P)  [13.] 

V  2  )  ftsincp 

where  the  desired  sidelobe  level  is  SL{ dB). 

Conventional  panoramic  surveillance  radars  are  normally  described  by  a  cosec- 
pattem  given  by 

g(cp)  =  S-— c  sc  ( cp )  [14.] 

The  angle  cp  is  the  elevation  angle  of  the  maximum  point  of  that  pattern.  A 
typical  value  of  cp  is  one-half  beamwidth.  The  directivity  factor  below  the  angle  (cpf;7f  + 
cp^/2)  is  represented  by  a  sine  function.  Figure  5  gives  examples  of  the  four  antenna  gain 
directivity  patterns  actually  used  in  CASE_ATTI  sensor  module. 

3.1.4  Target  Cross  Section 

The  mean  target  cross  section  a  depends  on  the  target  size  and  shape,  the  radar 
wavelength  and  the  target  aspect  angle.  Radar  targets  such  as  aircraft,  satellites,  missiles, 
and  ships  can  be  qualified  as  point  targets  when  viewed  by  conventional  radars  at 
appropriate  ranges.  A  table  can  be  used  to  describe  a  at  a  given  frequency  versus  0'  and 
(J)'  so  that  the  measured  a  can  be  used.  Otherwise,  the  model  used  to  describe  the  missile- 
type  targets  is  (Ref.  9)  given  by: 
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<j  =<rnjcos[\/<t>'2  +0'2  j  +  ^sin^Jty 


k'2+0'2 


[15.] 


where 

an:  nose-on  cross  section  (m2), 

elevation  angle  between  the  line  of  sight  and  the  target  speed  vector, 

0':  bearing  angle  between  the  line  of  sight  and  the  target  speed  vector, 

L:  length  of  the  target  (m), 

D:  diameter  of  the  target  (m). 

For  aircraft  and  ships  at  closer  ranges,  a  more  sophisticated  function  of  0'  and  4/ 
has  to  be  used.  A  target  cross  section  model  for  aircraft  is  given  in  the  SURSEM  model 
(Ref.  4). 


3.1.5  Atmospheric  Attenuation 

Depending  on  their  frequency,  electromagnetic  waves  at  RF  are  more  or  less 
affected  by  oxygen  and  water  vapor.  Part  of  their  energy  is  absorbed  and  another  part  is 
scattered  by  the  encountered  particles.  The  one-way  total  atmospheric  loss  (dB)  is  given 
by 


Lam  ~  Ko  +Rs(y  f,c+yr+  Jic  ) 


[16.] 


where 

L  :  losses  due  to  oxygen  and  water  vapor  absorption, 

Rs :  slant  range  to  the  target  (km), 

7f/c:  attenuation  coefficient  (dB/km)  caused  by  fog  or  clouds, 
7r:  attenuation  coefficient  (dB/km)  caused  by  rain  or  snow, 
7ic:  attenuation  coefficient  (dB/km)  caused  by  ice  clouds. 
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The  values  of  various  attenuation  coefficients  as  well  as  the  details  to  calculate  L* 0 
are  given  in  Ref.  2. 

3.2  System  Noise 

Natural  noise  is  the  ultimate  limiting  factor  in  the  detectability  of  radar  signals.  The 
origins  of  this  noise  lie  both  within  the  radar  and  in  the  radar  environment.  The  noise 
power  in  the  receiving  system  can  be  expressed  in  the  radar  equation  in  the  form  of  a 
system  noise  temperature  Ts  given  by  (Ref.  8,  p.  152) 

Ts  =  Ta  +  Tt(Lr  - 1)  +  LrT0{Nf  - 1)  [17.] 


where 

Tt :  thermal  temperature  of  the  transmission  line  taken  to  be  290  K, 

Lr:  receive  loss  (taken  to  be  1/2  Lt/r), 

Nf:  receiver  noise  figure  (RF  preamplifier,  mixer,  IF  amplifiers), 

Ta:  noise  power  received  by  an  antenna  from  natural  external  radiating  sources,  and 
To  =290  K. 

As  suggested  by  Blake,  (Ref.  8,  p.  172),  Ta  is  given  by: 

ra=(0.8767;'-254)/La+290  [18.] 

where  Ta'  is  obtained  from  Fig.  4.8  of  Ref.  8  and  La  is  the  dissipation  loss  within  the 
antenna.  For  reflector  antennas  and  active  arrays,  the  assumption  is  that  La  =  1  (Ref.  8,  p. 
385). 


The  thermal  noise  power  (watts)  at  the  output  of  the  IF  stage  is  defined  as 
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where  k  is  the  Boltzmann  constant  (1.3  8x1  O'23  J/K)  and  BN  is  the  noise  bandwidth  of  the 
receiver  in  hertz.  The  actual  noise  bandwidth  varies  with  the  transfer  function  of  the  filter 
or  set  of  filters  prior  to  detection  and  with  the  type  of  waveform.  When  pulse  compression 
is  used,  the  noise  bandwidth  is  taken  as  the  extended  bandwidth  of  the  pulse  compression 
filter  {Bn=  \!x'  where  %'  is  the  compressed  pulse  width). 


3.3  Signal-to-Noise  Ratio  (IF  Stage) 

The  signal-to-noise  ratio  at  the  output  of  the  IF  stage  is  given  by  the  ratio  of  [2] 
and  [19]  divided  by  some  additional  losses:  range-gate  loss  including  average  straddling 
loss  ( Lg ),  limiter  and  quantization  loss  (L/)  and  matching  loss  (Lm).  The  S/N0  is  given  by 


S  =  Pr 

N0  N'LX 


[20.] 


where  Lx  -  LgLlLm.  In  our  simulation,  we  assigned  the  following  typical  values  for  the 
losses: 


Lg  (dB)  =  2  (Ref.  8,  p.414); 

Li  (dB)  =  2  (Ref.  10,  p.  477),  and 


Lm(dB)  =  1  (Ref.  8,  p.  334). 
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FIGURE  5  -  Antenna  pattern  directivity 

a)  Omnidirectional 

b)  Gaussian:  ((pb=10°,  direction  20°) 

c)  Sine:  ((pb=10°,  direction  25°) 

d)  Cosec:  ((pb=20°,  direction  10°) 


-90 


P154344.PDF  [Page:  37  of  86] 

UNCLASSIFIED 

21 

4.0  INTERFERENCE:  CLUTTER  AND  JAMMING 

Clutter  and  jamming  are  limiting  factors  for  the  detection  of  targets.  In  this 
chapter,  the  characteristics  of  these  interferences  are  discussed  from  the  point  of  view  of 
the  radar  detection  range  performance;  mean  intensity  and  spectral  description  are  given. 
This  chapter  aims  at  computing  the  signal-to-clutter  and  signal-to-jamming  ratios  at  the 
output  of  the  IF  stage  just  before  signal  processing  takes  place  at  the  video  stage. 

4.1  Clutter 

Echoes  from  the  sea  and  precipitation  can  produce  very  significant  radar 
interfering  signals.  In  many  cases,  these  constitute  the  dominant  radar  performance 
limiting  factor.  As  ducting  allows  radiation  to  reach  targets  well  beyond  the  horizon,  it 
similarly  makes  possible  backscattering  from  beyond  the  horizon.  Because  ducting  effects 
on  target  echoes  are  important  to  describe  the  radar  long-range  detection  capability, 
ducting  must  also  be  considered  in  the  calculation  of  interference. 

The  difficulty  of  calculating  the  effects  of  ducting  on  clutter  reflectivity  is  that 
these  effects  have  just  been  recently  documented  and,  although  some  models  have  been 
proposed,  no  simple  model  appropriate  for  our  type  of  radar  model  has  been  satisfactorily 
validated.  Nevertheless,  we  contend  that  it  is  preferable  to  adopt  a  crude  representative 
model  than  to  simply  ignore  ducting  effects.  For  sea  clutter,  we  propose  an  empirical 
extension  based  on  clutter  data  reported  in  the  literature,  to  the  widely  accepted  Georgia 
Institute  Technology  (GIT)  model  (Ref.  11).  For  precipitation  (or  weather)  clutter,  we 
use  the  model  described  in  CARPET  (Ref.  2)  which  is  appropriate  for  a  surveillance  radar. 

4.1.1  Sea  Clutter  Power 

In  any  given  range  cell,  CR,  the  mean  power  received  from  sea  backscattering  (at 
the  input  of  the  RF  stage)  is  given  by 
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P  =  K. 


R: 


[21.] 


where 


Kr  is  the  radar  constant  as  defined  in  [3], 

g  is  the  antenna  directivity  factor  in  the  direction  of  the  illuminated  sea  patch  (unitless), 

Rsc  is  the  slant  range  of  the  illuminated  patch  (m), 

Go*  is  the  reflectivity  (also  known  as  normalized  clutter  RCS)  of  the  sea  surface  per  unit 
area  in  the  cell  (m2  /  m2), 

A  is  the  area  of  the  illuminated  sea  patch  (in  m2),  which  is  given  by:  A  =  RscQbCR, 

Gb  is  the  3-dB  horizontal  beamwidth  (expressed  in  radians),  and 
Cr  is  the  range  cell  (m). 


In  G0*,  the  asterix  means  that  the  reflectivity  incorporates  ducting  effects.  Under 
ducting,  large  sea  backscatter  energy  can  be  received  from  beyond  the  horizon  as  energy  is 
reflected  against  the  earth  over  the  horizon.  EREPS  (Ref.  12)  provides  a  simple  means  to 
compute  the  reflectivity  under  evaporation  ducting  but  the  model  is  known  to  significantly 
underestimate  clutter  reflectivity  under  most  conditions.  Here,  the  chosen  approach  is  to 
develop  a  new  empirical  model  based  on  published  data  (Refs.  13-15).  The  general 
expression  for  G0 *  when  expressed  in  dB  is 


G;(dB)  =  G0(R;c)  for  R'sc<Rtr 

r: 


:Go(/U-10Plog 


/? 


■3  (K-K)  f°r  Kc>K 


[22.] 


where 

Go  is  the  basic  reflectivity  as  given  in  textbooks  (dB), 

R*sc  is  the  terrestrial  range  of  the  clutter  patch  (km), 

Rtr  is  a  transition  range  from  where  the  empirical  extension  applies  (km), 
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P  is  the  reflectivity  decrease  coefficient,  and 
3  is  a  clutter  propagation  loss  (dB/km). 

For  G0,  the  GIT  model,  as  described  in  (Ref.  11),  is  used  since  it  is  believed  by 
several  authors  not  to  be  contaminated  by  refraction  effects.  The  model  is  valid  for 
frequencies  between  1  and  100  GHz.  The  transition  range  Rtc  is  taken  to  be  the  range  at 
which  the  grazing  angle  under  ducting  begins  to  differ  significantly  from  the  one  under 
standard  conditions.  Although  ray  tracing  could  provide  an  accurate  value  of  Rtc,  the  use 
of  a  rough  estimate,  applicable  for  a  shipbome-surface  sensor,  is  suggested  here  for 
convenience  considering  the  approximate  nature  of  the  model  extension.  When  there  is  no 
ducting  or  when  ducting  can  be  neglected,  (a  minimum  duct  height  of  1  m  is  necessary  to 
declare  a  ducting  condition)  Rtr  corresponds  to  the  horizon  range  between  the  antenna  and 
the  effective  wave  height;  this  is  the  normal  limit  of  applicability  of  the  GIT  model. 
Mathematically,  the  transition  range  can  thus  be  written  as 

Rtr  =  Rh  if  hduct  <  1  m 
Rtr  =  R,c  if  hduct>  lm 


with 


*»= 3. 57( 


[24.] 


and 


^rc  =  7- 1 — 3. 6  log  hduct 


[25.] 


where 

ha:  antenna  height  (m); 

he  =  0.375>va,  where  wh  is  the  significant  wave  height  (see  Ref.  10,  p.  273); 
hduct.  duct  height  (m). 
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Data  published  in  the  literature  indicate  that  P  can  vary  from  4  to  about  1  when  the 
conditions  go  from  standard  to  optimal  ducting  conditions.  The  optimal  duct  height  (i.e. 
for  which  surface  energy  is  maximum)  depends  on  the  radiated  frequency.  It  can  roughly 
be  estimated  by 


W)=37-8-69W) 


[26.] 


where/is  the  frequency  in  GHz.  Simple  linear  relationships  are  used  to  estimate  P  for  low 
sea  state,  Ssea  (sea  state  2  and  below),  based  on  the  data  found  in  the  literature.  For  higher 
sea  states,  a  correction  is  applied.  Consequently,  P  can  be  expressed  as 


with 


P  =  P° 
P  =  P° 


if  Ssea<  2 

~(Ssea-2)  ,  i 
4 


if  Ssea>  2 


[27.] 


3 

P  —  4 —  T  induct  if  induct  —  i^ducl(opt) 

duct(opt) 

3 

P  =  1  +  op.  _  ,  (fduct  ~  Kuct{opt) )  if  h duct  >  Kucllopt) 

nducl(opt) 


[28.] 


For  frequencies  less  than  10  GHz,  reported  data  show  a  slightly  greater  decrease  in 
reflectivity  with  range  than  that  given  by  the  second  term  of  equation  [27].  A  clutter 
propagation  loss,  S,  has  then  been  incorporated  to  match  the  published  data  better.  When 
ducting  is  negligible  or  absent,  3  is  used  to  describe  the  steep  decrease  of  reflectivity 
beyond  the  horizon.  The  parameter  3  is  given  by 
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3  =  3.6  for  hduct<  1m 

3  =  l±llog[_o.9(/)  +  9]  for  hduct>  1  m  and/ <10  GHz  [29.] 

3  =  0  for  >  1  m  and  /  >  10  GHz 


Figure  6  shows  a0*  versus  range  with  and  without  ducting  for  a  frequency  of  9.8 
GHz  using  vertical  polarization.  The  antenna  height  is  25  m  and  the  duct  height  is  10  m. 
Wind  speed  and  wave  height  are  4.5  m/s  and  1  m,  respectively. 


Sea  clutter 


FIGURE  6  -  Ducting  effects  on  sea-clutter  reflectivity 

4.1.2  Sea-Clutter  Power  Distribution  Skewness 

The  models  developed  for  the  generation  of  the  clutter  false  alarms  (see  Section 
7.2)  and  for  the  estimation  of  the  clutter  distribution  loss  ( Ld )  in  the  clutter  filtering 
process  (see  Section  5.2)  require  estimates  of  the  skewness  of  the  clutter  power 
distribution  with  respect  to  Rayleigh.  For  medium-to-high  resolution  radars,  distributions 
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of  clutter  return  have  a  longer  tail  than  Rayleigh  distributions  for  large-amplitude  clutters; 
the  median  can  differ  significantly  from  the  mean  value.  Unfortunately,  the  literature  on 
the  subject  is  insufficient  to  safely  establish  the  relationship  between  clutter  distribution 
and  the  depending  radar  and  environmental  parameters.  However,  representative  guesses 
can  be  made  based  on  measurements  recently  published. 

Let  us  define  the  product  beamwidth-pulselength,  X  =  Qb %'  (in  degree-jls),  to 
describe  radar  resolution.  According  to  the  published  data  (Ref.  10,  Fig.  7.14  and  Refs. 
16-17),  non-Rayleigh  distributions  are  obtained  for  X  <  4  °-j!s  when  the  pulselength  is 
shorter  than  4  |is.  For  greater  X  and  pulselength  values,  Rayleigh  distributions  are 
assumed. 

To  characterize  non-Rayleigh  distributions,  the  Weibull  formulation  is  used. 
Weibull  distributions,  which  are  described  with  only  a  few  parameters  (which  can  be 
related  to  the  physical  conditions  of  operation),  have  proven  to  be,  in  general,  better  than 
log-normal  distributions  to  describe  the  statistical  behavior  of  clutter  signals.  With 
Weibull,  the  mean  and  median  reflectivity  are  adapted  from  Ref.  18  as: 

ao°  =Oo-1.6a-l01°g(r(l+cO)  [30.] 

where  of ,  the  median  reflectivity,  and  Go  ,  the  mean  reflectivity,  are  both  expressed  in 
dB,  r  is  the  gamma  function  and  a  is  the  unitless  Weibull  parameter.  As  detailed  in 
Chapter  7,  the  parameter  a  characterizing  the  Weibull  distribution  can  be  used  in 
conjunction  with  the  calculated  mean  reflectivity  (previous  section)  to  estimate  the 
distribution  of  clutter  false  alarms.  In  addition,  the  median-mean  difference  is  used  to 
estimate  the  degradation  of  the  S/I  improvement  of  the  clutter  filtering  process  (called 
distribution  loss)  (see  Section  5.2).  The  parameter  a  depends  on  the  radar  characteristics 
(e.g.  resolution,  frequency,  polarization)  and  the  environmental  conditions  (wind  and 
grazing  angle).  For  medium-to-high-resolution  radars,  as  defined  above,  values  of  a  are 
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estimated  in  the  radar  domain  (R,Q)  and  stored  in  the  clutter  map  file  together  with  the 
calculated  reflectivities. 

From  experimental  results.  Ref.  19  shows  that  a  varies  roughly  linearly  with  the 
logarithm  of  wind  speed.  For  our  model  we  use 

aref  =  max^miH ,  3  log  ws)  [31 .] 

where  ws  is  the  wind  speed  in  m/s.  This  relationship  is  valid  at  X-band  (9.8  GHz),  upwind 
conditions  and  at  3°  grazing  angle.  To  obtain  the  a  value  in  any  conditions,  we  use  a 
relationship  between  a  and  the  magnitude  of  the  mean  backscatter.  Reference  19  shows 
that  a  varies  linearly  with  the  magnitude  of  the  mean  backscatter  when  expressed  in  dB. 
Consequently,  the  effective  a  parameter  can  be  given  by: 

a  =  max  (a^ , 0. 066a0d  —  aref )  [32.] 

where  God  is  the  difference  in  dB  between  the  evaluated  mean  reflectivity  g0  for  the 
conditions  under  study  and  the  Go  under  the  referenced  condition  stated  above.  The 
difference  God  must  then  be  precomputed  in  addition  to  Go*  and  stored  in  the  clutter  table. 

For  increasing  range  (or  decreasing  grazing  angle)  the  resolution  cell  widens, 
which  could  logically  contribute  toward  Rayleigh  distribution.  Reference  20  shows  how  a 
is  likely  to  decrease  with  increasing  resolution  cell  volume.  This  further  correction  of  a 
was  not  adopted  as  some  tests  made  using  various  measurements  found  in  the  literature 
showed  that  in  most  cases  the  resulting  a' s  are  within  the  margin  of  uncertainty  of  [31]  and 
[32].  In  our  model,  the  decrease  of  a  versus  range  is  then  essentially  governed  by  the 
decrease  of  G0  with  range  (i.e.  with  decreasing  grazing  angle). 

Some  authors  suggest  that  a  can  get  smaller  than  1  when  the  distribution  tends  to 
Rayleigh.  As  explained  in  Chapter  7,  for  the  generation  of  clutter-induced  false  alarms,  it 
is  convenient  to  have  this  number  adjustable  by  the  user.  While  a  small  variation  of  amin 
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can  affect  false  alarm  density,  it  has  no  significant  impact  on  the  clutter  filtering 
distribution  loss. 

The  method  described  above  for  estimating  the  skewness  of  the  sea-clutter 
distribution  is  based  on  clutter  statistics  which  apply  for  vertical  polarization  (VP) 
exclusively,  as  they  essentially  describe  the  Bragg  scatter  component  of  clutter  (Ref.  21). 
Unfortunately,  the  very  few  clutter  data  published  so  far  in  the  literature  for  horizontal 
polarization  (HP)  prevent  us  from  deriving  a  similar  model  for  HP  radar  emissions. 
Although  mean  backscatters  are  weaker  in  HP,  in  general,  distributions  have  longer  tails 
for  the  same  conditions.  Furthermore,  the  median-to-mean  difference,  a  skewness 
indicator,  has  often  been  found  to  be  maximum  at  sea  state  2  and  to  decrease  very  slightly 
from  sea  state  2  to  sea  state  5.  Consequently,  in  the  simulation,  until  more  information  is 
found  in  the  literature,  for  horizontal  polarization,  we  propose  to  use  a  high  value  of  a, 
a^,  ,  at  close  range  (reference  range  corresponding  to  a  grazing  angle  of  3°)  for  any 
wind  speeds  and  to  perform  a  smooth  linear  decrease  of  a  so  that  a  reaches  amin  at  the 
range  where  the  received  clutter  power  gets  down  to  2No  (i.e.  twice  the  thermal  noise 
power  (see  [19]);  a ^  is  arbitrarily  taken  to  be  3.0,  which  corresponds  to  a  relatively 
adverse  condition. 

4.1.3  Sea-Clutter  Spectrum 

The  sea-clutter  spectrum  width  and  shape  are  the  results  of  the  spread  of  velocities 
of  the  elemental  scatterers  of  the  sea  surface.  The  width  of  this  spectrum  is  an  important 
factor  in  assessing  the  performance  of  clutter  rejection  techniques.  A  Gaussian  spectrum 
shape  with  a  standard  deviation  as  yields 


W)  =  4exp 


2 


[33.] 
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The  mean  p,c  is  a  function  of  relative  motion  between  the  radar  platform  and  the 
clutter.  The  standard  deviation  of  the  clutter  spectrum  is  a  measure  of  the  bandwidth  of 
the  clutter  spectrum.  The  spectral  distribution  of  the  clutter  is  a  function  of  various  factors 
such  as  platform  motion,  internal  clutter  motion,  system  instabilities  and  transmitter  drifts. 

The  radar  platform  motion  creates  two  effects.  The  first  is  the  shifting  of  the  center 
frequency  of  the  clutter  spectrum  as  a  function  of  the  antenna  pointing  direction  relative  to 
the  platform  motion.  The  second  is  the  broadening  of  the  clutter  spectrum  resulting  from 
platform  motion  normal  to  the  antenna  pointing  direction.  The  finite  beamwidth  of  the 
antenna  illuminates  an  area  on  the  ground  or  sea,  which  contents  components  having  a 
different  Doppler  velocity  with  respect  to  the  radar.  The  Doppler  frequency  corresponding 
to  the  closing  velocity  between  the  antenna  and  the  center  of  the  clutter  patch  is  given  by 

/i=20J-jcos0asin\|frf  [34.] 


where 

Vj  is  the  ship  velocity  (m/s), 

0a  is  the  azimuth  angle  between  the  platform  velocity  vector  and  the  direction  of  the 
antenna  beam, 

'Prf  is  the  depression  angle  for  a  shipbome  antenna. 

There  can  also  be  a  Doppler  shift  due  to  the  average  radial  velocity  of  the  waves 
relative  to  the  radar,  which  is  typically  a  few  knots  (0.5  m/s).  The  mean  Doppler  shift  jlc  is 

Vc=fi-K*  P5J 


where  km,  is  roughly  1/A,.  Some  radars  compensate  for  the  platform  motion  and  in  this  case 
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The  width  of  the  clutter  spectrum  is  the  sum  of  several  components:  wind  shear, 
finite  antenna  beamwidth,  radar  system  instabilities,  etc.  The  scatterers  forming  the  clutter 
return  often  move  in  the  wind  causing  a  frequency  spreading  of  (Ref.  10) 


<*i  = 


OS(j)„ 

2% 


[36.] 


where  Gi  is  expressed  in  Hz  and 
ws  is  the  wind  speed  (m/s),  and 

(jV  is  the  direction  of  the  wind  (with  respect  to  true-north). 


If  the  radar  antenna  is  rotating,  an  additional  <Jr  resulting  from  the  antenna  rotation 
has  to  be  considered,  where  (Ref.  10,  eq.  9.76) 


<*r 


a 

5.356, 


[37.] 


where  a  is  the  antenna  rotation  rate  in  rd/s  and  0t  is  the  azimuth  beamwidth  (in  rd). 
Then,  the  resultant  standard  deviation  of  the  sea-clutter  spectrum  is  given  by 


°s  =Vaf  +0fr- 


[38.] 


4.1.4  Weather-Clutter  Power 

Raindrops  and  other  forms  of  precipitation  are  significant  reflectors  of  radar  waves 
and  thereby  can  produce  an  appreciable  clutter  echoes.  The  power  received  from  such 
echoes  at  the  input  of  the  receiver  is  given  by 


[39.] 
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where  Kr  is  as  previously  defined,  Rwc  is  the  slant  range  of  the  target  (m)  and  d  represents 
the  cross  section  of  the  volume  cell. 

For  a  surveillance  radar  where  the  beam  is  at  the  horizon,  the  model  described  in 
CARPET  (Ref.  2)  is  adequate.  The  cross  section,  ft,  is  given  by 

$  =  $0RcQbdltCR  [40.] 

with 

Rc  -  mod(/?,,  2^-)  [41.] 


and  where: 

t %  is  the  weather  reflectivity  (m2/m3), 

R,  is  the  terrestrial  distance  between  the  target  and  the  antenna  (m), 

Cr  is  the  range  gate  (m)  (=  ct72,  where  c  is  the  speed  of  light  and  x'  is  the  compressed 
pulsewidth), 

dh  is  the  height  of  the  volume  cell. 

The  weather  reflectivity  is  a  function  of  the  precipitation  rate,  pr  (mm/h),  and 
the  frequency,/ (GHz).  It  is  computed  using  the  general  expression: 

=  [42.] 

where 

K  —  1  x  lO^48  when  /<6GHz 
K  =  13X10-48  when  />  35  GHz 

For  6  GHz  </<  35  GHz,  K  is  obtained  from  a  linear  interpolation  between  7  x  10'48  and 
13  x  10'48.  This  formula  is  valid  for  rain-type  precipitation.  Snow  precipitation  rate  must 
then  be  converted  in  rain  equivalent  precipitation  rate.  For  dry  snow,  the  rain  equivalent 
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rate  is  accepted  to  be  10  times  less  than  the  actual  snow  precipitation  rate,  while  for  wet 
or  melted  snow  the  ratio  is  expected  to  be  less.  A  ratio  of  5  is  suggested  in  the  latter  case. 

The  height  of  the  cell  volume  is  given  by: 

dh  ~  min  (hcloud,  Rc$b )  [43.] 

where  hcioud  is  the  height  of  the  cloud  base  (height  of  the  rain  volume)  and  <\>b  is  the 
elevation  beamwidth. 

4.1.5  Weather-Clutter  Spectrum 

The  relative  motion  of  the  reflectors  produces  a  spread  of  velocities  and 
consequently  a  spread  of  Doppler  frequencies.  A  Gaussian  assumption  is  used  to 
approximate  the  weather  clutter-spectrum  shape.  Its  standard  deviation  (in  Hz)  is 
determined  by  (Ref.  10) 

^  _ _ _ 

=^VCTL +  <*L  t44-] 


with 


G shear  =  0. 42  X  4  X  10~3  i?  [45.] 

where 

(7  Shearis  the  wind  shear  standard  deviation  (Hz) 
a  lurb  is  the  standard  deviation  of  the  turbulence  (Hz) 

An  additional  deviation  due  to  antenna  rotation  is  added  as  in  the  case  of  clutter 


sea  spectrum. 
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4.1.6  Signal-to-CIutter  Ratio 

The  signal-to-sea-clutter  and  the  signal-to-weather-clutter  ratios  at  the  output  of 
the  IF  stage  are  given  by  the  ratio  of  P/Psc  and  P/Pm  both  affected  by  the  receiver  losses, 
Lc\ 


S  Pr 
c  ( Psc  +  PjLc 


[46.] 


where  Lc(dB)  =  Lx  + 1  dB,  i.e.  1  dB  more  than  the  loss,  Lx,  used  when  only  system  noise 
is  present,  as  slightly  greater  losses  are  expected  in  the  presence  of  clutter  (Refs.  8  and 
10). 

4.2  Jamming 

Only  standoff  noise  sidelobe  radar  jamming  is  considered  in  this  model.  In  standoff 
jamming,  the  jammer  stands  off  from  the  radar  at  a  fixed  position.  This  noise  jamming  has 
the  same  effect  on  the  target  detection  as  does  natural  noise;  the  only  essential  difference 
being  its  magnitude. 


The  power  of  the  jamming,  Pj  is  given  by 

P  -  Pj,GrGjLslBX2F f 
j  ( 4n)2  R]LrLatmBj 


[47.] 


where 

Pj, :  transmitted  jammer  power 
Gr:  radar  receiving  antenna  gain 
Gj\  jammer  antenna  gain 

Lsi.  azimuth  sidelobe  radar  antenna  level  with  respect  to  main  beam 
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B:  radar  beamwidth 
Ff.  one-way  propagation  factor 
Rf.  jammer  range 
Lr:  receive  loss 

Latm:  one-way  atmospheric  losses 
Bf.  jammer  bandwidth 


4.2.1  Signal-to-jamming  Ratio 

The  signal-to-jamming  ratio  at  the  output  of  the  IF  stage  is  given  by  the  ratio  of 
P/Pj  affected  by  receiver  losses  such  as 


l_P^J_ 

J  "  Pj  Lj 


[48.] 


where  Lj  =  Lx  as  in  the  case  of  system  noise. 
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5.0  RADAR  SIGNAL  PROCESSING 

The  performance  of  the  radar  detection  process  depends  on  the  signal  processing 
capability  to  reject  the  undesired  interference.  The  objective  of  this  chapter  is  to  compute 
the  S/I  ratio  at  the  output  of  the  video  stage  where  signal  processing  is  applied.  The  signal 
processing  yields  enhanced  signal  power  relative  to  clutter  and  noise  power.  This 
enhancement  is  termed  an  improvement  factor. 

5.1  Improvement  Factor  Calculation 

Velocity  filtering  for  modem  radar  systems  have  been  made  mainly  through 
moving  target  indicators  (MTI)  and  with  pulse  Doppler  radars.  MTI  is  normally  used  with 
low-pulse-repetition-frequency  (low-PRF)  radars.  The  MTI  is  usually  a  single  filter  that 
rejects  clutter  by  means  of  a  notch  in  its  passband  centered  on  the  clutter  Doppler 
spectrum.  Pulse  Doppler  is  usually  associated  with  high-PRF  radars.  A  contiguous  bank  of 
narrow-band  filters  is  used  to  detect  moving  targets  outside  of  the  clutter  spectrum.  This 
section  presents  the  improvement  factor  of  M-pulse  MTI  cancelers  as  well  as  the 
improvement  factor  obtained  with  a  coherent  integrator  (bank  of  narrow  filters). 

Both  the  MTI  and  the  coherent  integrator  can  be  modeled  as  a  transversal  filter 
where  the  output  y(t )  consists  of  a  weighted  sum  of  the  input  x(t)  such  as 

y(  0  =  X  wnx(J  ~  nT )  t49-] 

n 

where  T  is  the  time  between  samples,  n  is  the  filter  tap  number  and  while  the  weights  wn 
may  be  complex. 

The  transversal  filter  described  by  [49]  depicts  only  one  filter  and  one  output. 
However,  multiple  independent  outputs  are  possible  by  selecting  appropriate  weights.  In, 
particular,  the  FFT  weights  form  n  filters  with  n  outputs.  The  signal  gain  of  a  particular 
filter  is  found,  provided  that  the  Doppler  frequency,  fd,  is  such  that  a  signal  exists 
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somewhere  between  the  crossover  points  of  adjacent  filters.  The  gain  of  the  ith  filter  is 
derived  by  considering  a  target  Doppler  fd  that  can  occur  with  equal  probability  in  the 
region 


f,  ~  B/2  <fd<fi  +  B/2 

where/-  is  the  center  of  the  response  of  the  ith  filter  and  B  the  filter  bandwidth. 

The  probability  density  function  for  fd  associated  with  the  ith  filter  is 

P(/,)=-jr  for  f-B/2<f^f,+B/2 
-  0  elsewhere 

The  received  signal  waveform  can  be  represented  by 

r(t)  =  s(t)  exp(y27t//)  [52.] 

where  s(t)  takes  into  account  target's  characteristics,  transmitted  waveform,  radar  stability 
antenna  scanning  and  so  on,  and  fd  is  the  Doppler  frequency  corresponding  to  the  relative 
velocity  between  radar  and  target.  If  the  Doppler  shift  is  the  dominant  effect  then  the 
autocorrelation  of  eq.[52]  is 


[50.] 


[51.] 


\j/(t)  =  exp(j2nfdx) 


[53.] 


The  expected  value  of  \|/  is 

£{V(t)}  =  ]  atUWfiWM  =  exp02lt/t)S1”(-^>  [54.] 

J  TZBl 


and,  consequently,  the  signal  covariance  matrix  M.s  becomes 
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Ms(k,l)  =  cxp(-j2n{fi/fr)(k-l)) 


sin(7t  (B/fr)(k-l)) 
n{B/fr)(k-l ) 


[55.] 


It  is  convenient  to  assume  a  Gaussian  clutter  spectrum.  When  the  clutter  spectrum 
departs  from  the  Gaussian  assumption,  a  distribution  loss  is  taken  into  account  in  the 
overall  processing  gain  computation  (next  section).  Assuming  stationarity,  the  (k,l)th 
element  of  the  clutter  covariance  matrix,  Mc,  is  obtained  by  taking  the  Fourier  transform 
of  [55],  which  results  in 

Mc(k,  l)  =  exp(-2(7tGs  /  fr)2(k- 1)1  -  j'2k([Lc/ fr)(k  -  /))  [56.] 


The  improvement  factor,  IF,  of  any  arbitrary  transversal  filter  specified  by  a  weight 
vector  w  can  be  computed  using 


IF  = 


wHMcw 


[57.] 


where  wH  is  the  transposed  complex  conjugated  of  the  weight  vector  associated  with  the 
filter  of  interest. 


5.1.1  MTI  Filters 

The  improvement  factor  of  an  n-stage  MTI  can  be  computed  by  using  [57].  The 
matrix  Ms  is  generated  using  [55]  and  by  letting  =  0.5 fr,  i.e.the  center  of  the  filter  is  at 
one-half  the  pulse  repetition  frequency.  The  covariance  matrix,  Mc,  is  generated  using 
[56],  with  the  average  Doppler  clutter  ]lc  brought  back  to  zero.  The  column  vector  w  can 
be  given  by  the  binominal  coefficients  with  alternating  signs,  i.e. 


for  /  =  1,2, ...,«  + 1 


[58.] 
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5.1.2  Bank  Of  Narrow  Band  Filters 

Here  we  consider  the  more  general  approach  of  passing  the  data  through  a  bank  of 
narrow-band  filters.  Since  the  noise  is  wideband,  its  energy  is  divided  among  various 
filters.  However,  each  target  return  has  its  energy  concentrated  at  one  specific  Doppler 
frequency.  A  bank  of  narrow-band  filters  actually  reduces  the  average  value  of  the  noise 
by  eliminating  many  of  its  frequency  components.  We  consider  here  two  classes  of 
Doppler  filter  bank,  the  weighted  Discrete  Fourier  Transform  (DFT)  and  Finite  Impulse 
Response  (FIR)  filter  banks. 

The  DFT  filter  bank  comprises  a  set  of  filters  that  have  identical  characteristics, 
except  for  the  mainlobe  position.  Identical  filters  can  be  synthetized  using  the  FFT  (Fast 
Fourier  Transform).  FFT  filters  represent  a  restricted  class  of  FIR  filters  having  weighting 
coefficients  that  are  not  optimum  because  they  cannot  be  defined  independently  for  each 
filter.  The  FIR  technique  is  more  flexible  than  the  DFT  technique  for  the  design  of  the 
filter  transfer  function.  This  is  because  its  weights  can  be  changed  for  different 
environments  thus  changing  the  sidelobe  characteristics  of  the  filter.  If  necessary,  the  FIR 
technique  allows  for  filters  that  are  non-identical  in  order  to  concentrate  the  low  sidelobes 
in  the  frequency  range  of  the  clutter  interference. 

The  improvement  factor  for  coherent  integration  (bank  of  filters)  is  computed  the 
same  way  as  for  the  MTI  filters.  Equation  57  is  used  with  the  weight  vector  w  generated 
using 

for  n  =  1,2,...,  A.  [59.] 

where  k  is  the  filter  index  and  N  the  number  of  filters.  The  signal  covariance  matrix  M,  is 
generated  using  [55],  with  the  bandwidth  B  given  by  1  IN.  The  center  of  the  ith  filter,/-,  is 
given  by 


w 


kn 


exp 


-j2n 


N 
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for 


[60.] 


5.1.3  Cascaded  MTI  and  Coherent  Integrator 

In  this  case,  the  improvement  factor  is  calculated  still  using  [57].  The  MTI  filter  is 
characterized  by  the  number  M  of  pulses  processed  and  by  a  weight  vector  a.  The 
coherent  integrator  is  characterized  by  the  number  N  of  pulses  it  processes  and  by  its 
weight  vector  b.  The  number  M  is  always  less  than  N.  The  cascaded  combination  of  filters 
can  be  evaluated  as  outlined  for  the  MTI  filters  or  for  the  bank  of  narrow-band  filters.  The 
only  difference  is  that  the  weight  vector  w  is  generated  the  following  way: 


w  = 


«1 

0 

0 

...  o 

a2 

*h 

0 

...  o 

'■ 

0-2 

al 

•  •  •  • 

aM 

a2 

0 

0 

aM 

: 

a\ 

0  0a 

0  0  0 


M 


a , 


a. 


bi 


L  N-l  J(AT-l) 


M  Jmx(A'-I) 


[61.] 


5.1.4  Improvement  Factor  with  Staggered  PRF 

We  use  here  the  results  of  Ref.  22  where  the  effects  of  staggered  PRF  on  the  MTI 
improvement  factor  performance  is  investigated.  They  conjecture  that  the  minimum  clutter 
output  of  a  staggered  PRF  MTI  system  is  bounded  by  the  performance  of  two  equivalent 
constant  PRF  MTI  systems.  One  of  them  has  a  PRF  equal  to  the  lowest  PRF  of  the 
staggered  PRF  system,  while  the  other  has  a  PRF  equal  to  the  highest  PRF  of  the 
staggered  system.  In  our  computation  we  take  the  average  of  both. 
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5.1.5  Data  Windowing 

The  sidelobe  levels  of  the  Doppler  filters  can  be  reduced  by  using  a  data  window 
function.  This  corresponds  to  weighting  the  input  time  samples  x(n )  by  a  weight  vector 
a(n)  such  as 


w{n )  =  a(n)x(n )  [62.] 

where  a(n)  can  be  chosen  from  a  variety  of  window  functions  (Ref.  23)  for  a  desired  filter 
shape.  The  improvement  factor  is  then  calculated  using  [57]. 

5.2  Clutter  Distribution  Loss 

In  modem  radars,  short  pulselengths  (high  resolution)  are  often  priviledged  in 
order  to  minimize  clutter  returns  as  the  returned  mean  power  is  directly  proportional  to 
the  resolution  cell  size.  However,  as  a  drawback,  sea-clutter  spectrum  widen  and  the 
amplitude  distribution  shows  longer  tail  (for  strong  returns)  with  increasing  resolution, 
which  makes  sea  clutter  more  difficult  to  filter  out. 

For  medium-to-high  resolution  radars,  a  loss,  called  clutter  distribution  loss,  is 
heuristically  estimated  to  account  for  the  degradation  of  the  improvement  factor,  IF,  in 
non-Rayleigh  clutter  when  sea  clutter  is  found  to  be  the  dominant  clutter  source.  The 
conditions  for  identifying  medium-to-high  resolution  radars  are  as  defined  in  Subsection 
4.1.2.  In  our  model,  we  assume  that  the  losses  presented  by  Schleher  (Ref.  24),  as  a 
function  of  IF  and  the  clutter  mean- to -median  ratio,  describe  the  radar  filtering  efficiency 
degradation. 

From  Fig.  1  of  Ref.  24,  we  come  up  with  the  approximate  expression  for  the 
clutter  distribution  loss  L/. 

Ld  =  Ld0  for  0  <  IF  <  20  dB 

=  Ld o  exp(-(/F  -20 )2  /  lOo)  for  I F  >  20  dB 


[63.] 
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This  relationship  is  believed  to  provide  realistic  values  of  Ld  for  Pd  vaying  from  0.5  to  0.8, 
the  normal  threshold  of  detection.  Lj0  can  be  approximated  by: 

Ld0  =  1 5(  1  -  exp(—  dmm  /8. 45))  -  5  [64.] 

where  d,Jim  is  the  difference  between  the  median  and  the  mean  clutter  return.  For  any  given 
radar  cell,  dmm  can  be  calculated  from  [30]  (Subsection  4.1.2)  using  the  Weibull  parameter 
a  and  the  mean  clutter  reflectivity  given  in  the  clutter  map  table. 

5.3  Processing  Gain 


The  processing  gain,  Gp,  of  the  video  stage  is  consequently  given  by 


[65.] 


An  example  of  computation  of  IF  is  given  in  Fig.  7  where  five  HR  filters  with  32  Hanning 
filter  taps  have  been  used.  The  IF  is  plotted  versus  Doppler  frequency  normalized  with 
respect  to  PRF.  The  oscillations  on  the  curve  are  caused  by  a  phenomenon  called  the 
scalloping  or  picket-fence  effect.  The  scalloping  effect  is  the  observed  dip  in  the 
improvement  factor  for  a  Doppler  frequency  lying  midway  between  two  filters. 


improvement  Factor  (t») 
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FIGURE  7  -  Improvement  factor  versus  normalized  Doppler  frequency 
for  a  bank  of  5  FIR  filters 
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6.0  TARGET  CONTACTS 

A  target  is  detected  when  the  received  signal  exceeds  a  pre-set  threshold  which 
depends  on  a  desired  probability  of  false  alarm  Pfa-  The  probability  that  the  threshold  be 
exceeded  due  to  the  presence  of  a  target,  i.e.  Pd,  depends  on  noise,  clutter  and  jammer 
statistics,  threshold  level  and  integration  of  pulses  as  the  beam  sweeps  the  target.  In  the 
previous  chapters,  we  calculated  the  necessary  components  to  evaluate  the  signal-to- 
interference  energy  ratio  resulting  from  the  reception  and  processing.  In  this  chapter,  the 
S/I  ratio  is  used  to  calculate  the  probability  of  detection  for  a  given  false-alarm  probability. 
Then,  “realistic”  target  contacts  measurements  are  generated  according  to  Pd  and  the 
accuracy  of  the  sensors. 

6.1  Calculation  of  the  Signal-to-Interference  Ratio 


The  equivalent  5/1  ratio  per  pulse  is  calculated  according  to  the  following 
expressions.  When  thermal  noise  is  dominant  (i.e.  N0'  >  2C  >  J  where  C  results  from  both 
sea  and  weather  clutter),  S/I  is  given  by 


5=_S_g  1 

I  No  P  Lcfar 


[66.] 


where  Lcfar  is  the  loss  associated  to  set  up  the  threshold.  When  thermal  noise  is 
dominant,  Ld=  1.  In  CFAR  system,  the  threshold  is  constinuously  adjusted  at  some  factor 
times  the  noise  estimated  from  a  second  reference  channel.  When  the  number  of  samples 
of  noise  in  the  reference  channel  is  reduced,  the  sampled  reference  noise  will  fluctuate 
appreciably  and  a  higher  threshold  must  be  set  up.  LCfar  corresponds  to  a  loss  in  signal 
detectability.  A  value  of  1.5  dB,  which  is  a  likely  average  value  for  operational  systems 
(Ref.  8),  has  been  assigned  to  Lcfar  ■ 


When  the  clutter  is  dominant  (i.e.  C>  N0'/2 ),  the  S/I  is  given  by 
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S_ 

I 


^CFAR^cd 


[67.] 


where  Lcd  describes  a  decrease  in  detection  sensibility  when  the  clutter  amplitude 
distribution  is  non-Rayleigh.  A  value  of  2  dB  has  been  assigned  to  Lcd,  which  correspond 
to  a  likely  trade-off  for  theshold  setting  in  order  to  optimize  detectability  while 
maximazing  clutter  rejection  (Ref.  18).  For  the  beam  shape  loss  LB  ,  the  conventional 
value  of  1 .6  dB  from  Blake  (Ref.  8)  is  applied. 


When  jammer  is  dominant  ( J  >  2(C+  N0')),  it  is  similar  to  thermal  noise. 


6.2  Calculation  of  the  Probability  of  False  Alarms 

In  radar  specifications,  the  mean  time  between  false  alarms,  tfa,  is  the  parameter 
that  describes  the  desired  false  alarm  rate  instead  of  the  probability  of  false  alarm,  which  is 
used  in  the  calculation  of  the  probability  of  detection.  As  discussed  in  Ref.  10,  the 
probability  of  false  alarm  Pfa  is  given  by: 


P  _  °-69 


[68.] 


where  n'  is  the  number  of  independent  opportunities  for  a  false  alarm  in  the  mean  time 
between  false  alarms,  tfa.  Defining  n0  as  the  number  of  opportunities  during  the  beam 
sweep  period  over  any  potential  point  target,  and  tmv  the  processing  time  required  to 
determining  whether  or  not  a  target  is  present,  n'  can  be  expressed  as: 


n'  =  n0  tfa  Jt^p 


[69.] 


where 


"  2n 


[70.] 
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and 


n0  -  NchNrc 


[71.] 


where 

tscan  is  the  time  of  one  radar  revolution  (s), 

Qb  is  the  horizontal  3  dB  beamwidth  (rad), 

Nch  is  the  number  of  radar  decision  channels  (unitless  integer),  and 
Nrc  is  the  number  of  range  cells. 


The  number  of  range  cells  is  given  by 


N 


-2 


[72.] 


where  fr  is  the  pulse  repetition  frequency  (Hz)  and  x'  is  the  compressed  pulse  length  (s). 
In  our  model,  the  number  of  decision  channels  per  cell,  Nch  is  taken  to  be  2  for  any  radar. 

6.3  Effective  Number  of  Pulses  Available  for  Noncoherent  Integration 

When  system  noise  or  jamming  are  the  dominant  sources  of  interference,  we 
assume  that  the  interfering  signal  is  Gaussian  and  uncorrelated.  We  also  assume  that,  in 
this  case,  clutter  filtering  stages  are  switched  off  by  the  operator.  Consequently,  the 
number  of  pulses  available  for  non-coherent  integration  is  the  number  of  pulses  that  hit  the 
point  target  during  a  beam  sweep,  which  is  given  by 


N n  -tswp  ■  fr  [73.] 

When  clutters  dominate  as  sources  of  interference,  the  interfering  signals  are 
correlated;  this  contributes  to  decrease  the  effectiveness  of  the  non-coherent  integration  of 
pulses.  To  account  for  the  integrator  performance  degradation,  we  define  an  effective 
number  of  pulses  available  for  the  integration,  Ne  ,  which  is  a  number  between  1  and  Nn 
depending  on  the  radar  and  clutter  characteristics.  As  explained  by  Nathanson  (Ref.  10 
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pp.  89-92),  the  nominal  effective  number  of  pulses,  here  defined  as  Nn' ,  can  be  understood 
as  the  number  of  independent  (or  uncorrelated)  pulses;  it  can  be  approximated  by 

W,'=[W,]  if  (or;/V>  0,2 
W„'=[l  +  5a;etM/J,  eke. 

where  X  is  the  wavelength  (m)  and  a/  is  the  standard  deviation  of  the  clutter  velocity 
spectrum  (m/sec). 

In  most  modern  radars,  the  sampling  rate  is  such  that  clutter  samples  are  highly 
correlated  and  thence  Nn'  is  very  small  (i.e.  1  or  2)  under  most  conditions.  However, 
frequency  agility,  which  is  often  employed,  is  a  means  to  decorrelate  clutter  samples. 
Therefore,  the  effective  number  of  pulses  available  for  non-coherent  integration  before  (or 
without)  MTI  processing,  hereafter  referred  to  as  N/  can  be  generalized  by  the  expression: 

N’  =  ma(N,/Nf,K)  [75.] 

where  Nfc  is  the  number  of  pulses  before  a  frequency  change  (as  given  in  the  radar 
specifications). 

MTI  filtering  contributes  to  further  correlate  incoming  interfering  signal  and  thus 
to  decrease  the  number  of  effective  pulses.  Assuming  that  the  N/  pulses  as  obtained 
above  are  as  uncorrelated  as  system  noise  (so  that  they  can  be  considered  like  noise 
pulses),  the  number  of  effective  pulses  at  the  ouput  of  the  MTI  filtering  can  be 
approximated  using  a  fit  on  theoritical  results  published  by  Trunk  (Ref.  25): 


N=- 


K 


1  +  0.53  Npc 


[76.] 


where  Npc  is  the  order  of  the  MTI  filter  (A^-pulse  canceler). 


P154344.PDF  [Page:  63  of  86] 


UNCLASSIFIED 

47 


Finally,  when  filter  bank  or  coherent  integration  is  used  instead  of  MTI  the 
effective  number  of  pulses  available  for  post-integration  is  (Refs.  8,  10): 


Ne  =  NJNc 


[77.] 


where  Nc  is  the  number  of  pulse  coherently  integrated. 


6.4  Calculation  of  the  Probability  of  Detection 

A  convenient  and  conceptually  simple  approach  based  on  the  widely  accepted 
Blake  procedure  (Ref.  8)  is  used  to  calculate  the  probability  of  detection.  In  Blake’s 
approach,  a  minimum  signal-to-noise  ratio  is  determined  to  achieve  a  given  probability  of 
detection  and  false  alarm.  Our  approach  differs  in  that  the  probability  of  detection  is 
determined  from  an  equivalent  ( S/I)  ratio  per  pulse. 

A  fairly  general  approximation  of  the  equivalent  S/I  per  pulse,  D0,  for  the 
detectability  of  steady  state  targets  is 


where 


1+  1  + 


Mjl 


X0  ={gfa  +  Sd)2 
gfa  =  2.367=4^—1.02 


1.23  R 


[78.] 


[79.] 

[30.] 

[81.] 


f  =  0.9(2Rd-l) 


[82.] 
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with 

Hn :  equivalent  number  of  pulses  integrated, 

£:  asymptotic  efficiency  of  the  envelope  detector, 

Ne :  number  of  integrated  pulses  (calculated  in  the  previous  section). 

This  formula  is  valid  over  the  range  0.1  <  Pd  <  0.9  and  for  10-12  <  Pfa  <  10”4  with 

a  maximum  error  of  0.7  dB  when  compared  with  theoretical  curves  (Figs.  2.4-2. 8  of  Ref. 
8).  The  equivalent  number  of  pulses  integrated,  Hn,  for  several  types  of  integrator  and 
two  signal  variation  models  is  given  in  Table  I  (Ref.  8,  p.  68). 

TABLE I 

Equivalent  number  of  pulses  integrated  for  different  kind  of  integrators 


SIGNAL  VARIATION 

INTEGRATOR 

Hn 

Gaussian 

Non-uniform  weight 

0.626  Ne 

Gaussian 

Uniform  weight 

0.558  Ne 

Gaussian 

Single-loop  feedback 

^ 1 

Gaussian 

Double -loop  feedback 

0.586  Ne 

Constant 

Non-uniform  weight 

Ne 

Constant 

Uniform  weight 

Ne 

Constant 

Single-loop  feedback 

1.4299  Ne 

Constant 

Double-loop  feedback 

1.3347  Neexp(-3/Ne) 

Table  II  gives  the  asymptotic  efficiency  £  for  the  three  envelope  detectors. 
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TABLE  II 

The  asymptotic  efficiency  of  three  envelope  detectors 


ENVELOPE  DETECTORS 

c 

Square  law 

1.0 

Linear  law 

0.915 

Logarithmic  law 

0.608 

The  probability  of  detection  Pd  is  obtained  by  an  iterative  procedure  to  an  accuracy 
of  1%  of  the  input  S/I  ratio.  Equation  78  is  valid  for  the  non-fluctuating  target  case.  For 
fluctuating  target  cases,  a  fluctuation  loss  in  dB  is  subtracted  from  D0  in  the  iterative 
process.  The  general  fluctuation  loss  for  the  five  Swerling  cases  that  are  normally 
considered  in  radar  systems  is  given  by  (Ref.  8,  p.75). 

£,=[(- to  Pj(l+W*JF  [83.] 

where  K  is  given  by  Table  III.  The  Swerling  case  3  is  used  when  a  missile-type  target  is 
considered.  For  large  targets,  like  aircraft,  the  user  can  select  either  cases  although  a  case 
1  is  recommended. 


TABLE  m 

Parameter  K  of  the  chi-square  family  to  describe  the  Swerling  cases 


TARGET  MODEL 

K 

Non-fluctuating 

— >oo 

Swerling  I 

1 

Swerling  II 

Ne 

Swerling  HI 

2 

Swerling  IV 

2  Ne 
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6.5  Generation  of  Measurements 

The  generation  of  realistic  "corrupted"  measurements  and  the  provision  of  the 
expected  measurement  error  (to  the  tracking  and  data  association  algorithms)  are  needed 
for  the  data  fusion  analysis.  Surveillance  radars  detect  targets  in  space  and  report  it  in 
(i?,9)  plane.  The  measurement  errors  indicate  the  accuracy  of  the  radar  output,  while  the 
resolution  indicates  the  capability  to  resolve  close  targets.  The  measurement  errors  in 
range  and  azimuth  depend  on  the  quality  of  the  received  signal  and  thereby  is  a  function  of 
the  target  location.  Simple  relations  between  the  measurement  errors  and  S/I  are  provided 
in  textbooks  and  are  generally  used  in  tracking  studies. 


The  mean-range  measurement  error  er  can  be  expressed  as 


[84.] 


while  for  the  mean  angle  measurement  error  £e  we  have 


[85.] 


where  Cr  is  the  range  cell  and  ( S/I)out  is  the  signal-to-interference  ratio  at  the  output  of  the 
threshold  stage  (Fig.2)  (i.e.  which  includes  the  effect  of  the  non-coherent  integration). 
Usually,  radar  specifications  give  Er  and  e9  for  a  specific  S/I  ratio  which  allows  to 
determine  constants  Hr  and  He.  When  they  are  not  supplied,  Hr  -  Hq.  =4  is  used  as  a 
default  value.  With  this  value,  the  mean  error  corresponds  to  the  resolution  when  ( S/I)ou , 
is  12  dB  (a  likely  nominal  detection  threshold). 

The  measured  target  position  is  generated  by  adding  random  errors  to  the  true- 
target  position  according  to 


Rmeas  =  Kue  +  N(°>  /(£«)) 


[86.] 
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=0r™  +  N(O,/(£e))  [87.] 

where  N  is  a  Gaussian  distribution  with  zero  mean  and  with  the  variances  taken  to  be 

/(£«)  =  (£*/  3)2  [88.] 

/(ee)  =  (£e/3)2  [89.] 

respectively,  so  that  99%  of  cases  (plus-or-minus  3  times  the  standard-deviation)  are 
within  plus-or-minus  the  average  value. 
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7.0  GENERATION  OF  FALSE  ALARMS 

For  the  study  of  sensor  data  correlation  and  fusion,  generation  of  realistically 
spread  out  false  alarms  is  as  important  as  the  predictions  of  detected  target  positions  since 
the  performance  of  fusion  algorithms  are  to  be  assessed  with  respect  to  their  capability  of 
establishing  target  tracks  out  of  false  alarms. 

Modern  radar  receivers  include  constant  false  alarm  rate  (CFAR)  devices  to  adjust 
the  detection  threshold  level  so  that  Pfa  remains  constant  in  the  presence  of  noise  or 
interference.  The  need  of  CFAR  arises  when  the  noise  level  in  the  receiver  varies  with 
time  in  some  unpredictable  manner,  when  clutters  are  significant  or  in  the  presence  of 
jamming.  In  the  present  radar  model,  we  attempted  to  describe  in  a  realistic  manner  the 
false  alarm  field  caused  by  the  system  noise  and  sea  clutter  in  the  instrumental  range- 
azimuth  domain  (R,0)  of  the  radar. 

In  the  calculation  of  target  Pd  ,  Pfa  is  assumed  to  be  kept  constant  in  any 
circumstances  at  the  cost  of  minor  S/I  losses,  called  CFAR  losses.  With  the  assumption  of 
a  constant  false  alarm  rate,  an  increase  of  noise-like  interferences  results  in  a  loss  of 
detectability  as  the  S/I  decreases  and  the  threshold  increases. 

False  alarms  caused  by  system  noise  occur  at  a  rate  corresponding  to  the  mean 
time  between  false  alarm  tfa  which  is  directly  related  to  Pfa  (see  Section  6.2).  Section  7.1 
describes  the  technique  used  to  generate  noise-induced  false  alarms. 

Sea  clutter  can  also  be  an  important  source  of  false  alarms  since  CFAR  systems 
can  hardly  adapt  to  the  clutter  statistical  behavior  which  may  substantially  differ  from  that 
of  Rayleigh  noise.  The  density  of  sea-clutter  false  alarms  depends  on  the  radar 
characteristics  (especially  resolution,  frequency,  polarization),  the  environmental 
conditions  (wind  speed)  and  the  grazing  angle.  The  number  of  sea-clutter  false  alarms  is 
normally  much  larger  than  the  one  due  to  noise.  These  false  alarms  are  not  uniformly 
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distributed  in  the  ( R,Q )  domain.  Their  occurrence  is  expected  to  decrease  with  range  (with 
decreasing  grazing  angle)  and  to  increase  in  upwind  direction. 

The  relationship  between  sea  clutter  false  alarm  probability,  Pfac  and  the  dependent 
factors  is  not  sufficiently  known  to  accurately  predict  Pfac  versus  range  and  azimuth  for  a 
diversity  of  radar  characteristics  and  sea  conditions.  The  literature  essentially  presents 
measurements  of  backscatter  coefficient  distribution  for  a  limited  set  of  environmental 
conditions  and  radar  characteristics,  and,  furthermore,  very  little  can  be  found  on  the 
capability  radars  to  filter  strong  clutter  returns.  In  spite  of  the  incomplete  knowledge  on 
sea-clutter-induced  false  alarms,  we  deem  that  it  is  preferable  to  use  a  heuristic  model  that 
provides  clutter-like  spreading  of  false  alarms  rather  than  using  a  uniform  spreading.  The 
heuristic  model  of  false  alarm  generation  uses  an  approximative  global  false  alarm  rate  that 
can  be  adjusted  by  the  user  and  it  features  a  spreading  scheme  which  is  representative  of 
published  experimental  results. 


7.1  Noise-Induced  False  Alarms 

False  alarms  can  occur  in  any  of  the  Nceus  resolution  cells  of  the  radar  in  the 
instrumental  range-azimuth  domain  (/?,„«, 0),  where 


M  -360 
c'"s  Q„cr 


[90.] 


•  A  detection  is  declared  in  a  cell  if  the  output  of  the  signal  processor  for  this  cell 
exceeds  a  threshold.  If  the  sensor  points  at  a  region  where  there  is  no  target,  then, 
detections  can  still  occur  due  to  noise  with  a  probability  given  by  Pfa-  Assuming  that 
the  events  in  the  cells  are  independent,  the  probability  of  having  m  false  alarms  in  the 
(Rinstfi)  domain  is  approximated,  using  (r\=NceusPfa),  by 


•  a  Gaussian  distribution  with  mean  and  variance  equal  to  rj,  when  r|  >  9; 
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•  a  uniform  distribution  when  r|  < 


0.69 

IN. 


ch 


t/(0,l)<Tl 

otherwise 


•  otherwise,  a  Poisson  distribution  with  mean  rj  is  used; 

(h)m  exp(-ri) 


P(nfo  =m)  =  - 


ml 


[91.] 


[92.] 


[93.] 


When  generating  the  Poisson  distribution,  the  "innovation  filter"  method,  as 
described  in  Ref.  26,  is  used  to  determine  the  number  of  false  alarms  for  each  scan:  nfa.  In 
this  technique,  the  cumulative  distribution  (either  Poisson  or  normal)  is  first  evaluated. 
Then,  a  random  number  from  0  to  1  is  drawn  from  a  uniform  distribution.  The  output  nfa 
is  the  number  of  event  for  which  the  cumulative  probability  corresponds  to  the  drawn 
number.  The  last  step  consists  in  determining  randomly  nfa  positions  in  the  (R,9)  domain 
using  an  uniform  generator. 

7.2  Sea-Clutter-Induced  False  Alarms 

For  sensor  data  fusion  study,  it  is  required  to  have  a  model  that  gives  a  realistic 
representation  of  the  spreading  of  sea-clutter-induced  false  alarms  as  opposed  to  uniformly 
distributed  false  alarms.  As  our  current  knowledge  on  sea-clutter  false  alarms  is  not 
sufficient  to  achieve  an  accurate  generation  of  false  alarms  as  a  function  of  the  dependent 
factors  (e.g.  radar  characteristics,  environment),  we  essentially  searched  for  a 
representative  spatial  variation  of  the  sea-clutter-induced  false  alarm,  Pfac,  using  likely 
values  of  global  probabilities. 

In  the  case  of  medium-to-high-resolution  radars,  sea-clutter  amplitude  distributions 
differ  from  the  Rayleigh  distribution;  this  is  mainly  due  to  the  significant  occurrence  of 
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very  large  clutter  returns,  substantially  larger  than  the  mean  backscatter  coefficients  <j0  as 
calculated  in  Subsection  4.1.1.  These  large  clutter  returns,  commonly  referred  to  as 
spikes,  are  responsible  for  the  formation  of  a  longer  tail  in  the  backscatter  distribution 
compared  with  Rayleigh  distribution.  Furthermore,  spikes  have  been  shown  to  have 
shifted  and  have  larger  spectrum  than  clutters  of  mean  magnitudes  (Ref.  16),  which  makes 
them  even  more  difficult  to  be  filtered  out  by  radars. 

Our  model  is  based  on  the  assumption  that  the  CFAR  devices  are  designed  to 
adjust  the  threshold  to  obtain  the  required  input  Pfa  assuming  that  clutters  are  Rayleigh 
distributed.  Then,  estimating  the  actual  clutter  distribution  for  a  cell,  one  can  compute  the 
probability  that  clutters  exceed  the  threshold  for  this  cell.  It  is  assumed  that  all  returns 
greater  than  this  threshold  (hereafter  called  spikes)  have  the  required  magnitude  and 
spectral  content  not  to  be  filtered  out  by  the  radar  clutter  rejection  stages. 

The  occurrence  of  spikes  has  been  shown  to  increase  with  radar  resolution  and  is 
thought  to  be  negligible  for  low  resolution  radars.  The  data  published  in  Ref.  10  (Fig. 
7.14)  and  Refs.  16-17  suggest  that  non-Rayleigh  distributions  can  be  obtained  for  pulse 
lengths  less  than  a  few  |is.  Defining  the  product  beamwidth-pulselength,  X,  as  in 
Subsection  4.1.2,  we  can  assume  that  spike-false  alarms  occur  for  X  <  4°-ps  when  the 
pulselength  is  shorter  than  4  ps.  Significant  occurrences  of  spikes  were  observed  by  Chan 
(Ref.  16)  in  vertical  polarization  and  Sekine  et  al.  (Ref.  17)  in  horizontal  polarization  for  X 
=  1  and  3.6°  ps,  respectively.  For  greater  X  and  pulselength  values,  clutter  distribution 
becomes  very  close  to  Rayleigh's  and  so  clutter-false  alarms  can  be  considered  negligible. 

Although  it  is  now  widely  recognized  that  clutters  are  better  characterized  by  a  Re¬ 
distribution  (Ref.  16),  Weibull  distributions,  being  characterized  by  fewer  parameters,  are 
more  convenient  and  they  offer  satisfactory  performances  (Refs.  16-18).  In  addition,  the 
Weibull  parameters  can  be  related  to  physical  conditions  of  operation.  Using  Weibull,  the 
probability  of  clutter  false  alarms  is  the  probability  of  having  a  return  exceeding  a 
threshold  Y : 
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Pfac  PweibullO',a') 


— ln(2) 


V^O median  J 


[94.] 


where  a  is  the  Weibull  parameter  and  <50median  is  the  median  backscatter  coefficient.  For  a 

=1  the  distribution  reduces  to  Rayleigh's.  In  the  calculation,  the  mean  a0  (Subsection 
4.1.1)  is  used  in  [94]  instead  of  o0median. 


Thence  for  any  given  sector  of  the  (R,Q)  domain,  the  probability  of  sea-clutter  false 
alarm,  Pfac,  is  calculated  with  [94]  using  a0  (mean)  and  a  as  given  in  the  sea  clutter  table 
and  the  threshold  F#? 

®  0  median  +  amin(101°g(ln(1///a))  +  1-6)  [95.] 

which  corresponds  to  the  required  threshold  to  obtain  Pjac=  Pfa  in  the  case  of  Rayleigh 
clutter  (amin  =  1).  Here  again,  the  mean  value  Ob  is  used  instead  of  the  median  value, 
a0 median-  The  probability  Pfac  can  then  be  evaluated  for  all  sectors  as  defined  in  the  clutter 
map.  The  technique  described  in  Section  7. 1  is  then  used  to  uniformly  generate  the  false 
alarms  in  sectors  from  their  respective  Pfac. 


Some  authors  (Ref.  15),  show  that  the  Weibull  parameter  a  may  be  smaller  than  1 
(e.g.  0.8)  at  long  ranges  or  for  low  mean  backscatters.  Consequently,  we  propose  to  let 
Pfac  user-modifiable.  At  the  same  time,  that  gives  the  user  a  control  of  the  global  false 
alarm  rate  density.  The  parameter  amin,  is  set  to  1  and  a  new  parameter  Kciutter  is  defined  so 
that 


p  _  PweibullO' j  Q) 
*fac  ' 


K, 


clutter 


[96.] 


where  KciutterC an  be  greater  or  equal  to  1. 
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For  sectors  where  a  gets  down  to  amin  or  where  the  received  sea-clutter  power  is 
of  the  order  of  the  system  noise  (2N0'  is  arbitrary  used),  one  can  assume  that  sea  clutter 
does  behave  like  systems  noise  and  so  Pfa  -  0  since  the  CFAR  is  assumed  fully  efficient  in 
this  case. 

The  false  alarm  generation  technique  hereabove  described  is  based  on  clutter 
statistics  which  apply  for  vertical  polarization  (VP)  exclusively,  as  they  essentially  describe 
the  Bragg  scatter  component  of  clutters  (Ref.  21).  Unfortunately,  the  very  few  clutter 
data  published  so  far  in  the  literature  in  horizontal  polarization  (HP)  prevent  us  to  derive  a 
similar  false  alarm  generation  technique  for  HP  radar  emissions.  Although  mean 
backscatters  are  weaker  in  HP,  in  general,  distributions  have  longer  tails,  and,  as  shown  by 
Chan  (Ref.  16),  the  spike's  spectrum  in  HP  is  shifted  more,  which  suggests  that,  for  a 
given  threshold  adjusted  above  the  mean/median  value,  much  greater  Pjac  would  be 
obtained  in  HP  than  in  VP. 

Calculations  performed  using  distributions  published  by  Olin  (Ref.  27)  indicated 
that,  in  HP,  the  probability  Pfac  can  be  several  orders  of  magnitude  greater  at  sea  state  2 
and  a  bit  less  than  10  times  greater  at  sea  state  5  as  compared  with  VP;  these  results  agree 
with  measurements  published  by  Ewell  (Ref.  28).  Furthermore,  in  contrast  with  VP,  in 
HP,  Pfac  is  expected  to  slightly  decrease  with  wind  speed  (and  probably  so  with  increasing 
mean  backscatters);  a  decrease  of  Pfac  of  less  than  an  order  of  magnitude  can  be  predicted 
from  sea  state  2  to  sea  state  5  based  on  measurements  published  by  Olin  (Ref.  27). 

Thence,  in  HP,  the  density  of  false  alarms  is  not  likely  to  decrease  with  range  (and 
with  the  magnitude  of  mean  backscatter)  as  fast  as  in  VP.  Furthermore,  the  global  density 
of  false  alarms,  which  does  not  vary  much  with  sea  state,  is  expected  to  be  significantly 
greater  than  in  VP  (likely  to  correspond  about  to  the  Pfac  at  very  high  sea  state  in  VP  in 
most  cases).  Consequently,  in  the  simulation,  until  more  information  is  found  in  the 
literature,  for  HP,  we  propose  to  use  a  relatively  high  maximum  PfaCmax,  at  close  range 
(reference  range  corresponding  to  an  angle  of  3°)  for  any  wind  speeds  and  to  perform  a 
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smooth  decrease  of  Pfac  up  to  the  range  where  the  parameter  a  reaches  amin,  or  where  the 
received  mean  clutter  power  reaches  2 No'.  At  further  ranges,  Pfac  is  taken  to  be  0.  For 
Pfacmax,  an  arbitrary  number  between  10'2  and  10'3  can  be  selected  by  the  simulation  user  or 
one  can  use  the  probability  corresponding  to  a  =  3.0  in  VP. 

The  two  histograms  of  Fig.  8  (a)  and  (b)  show  a  comparison  of  the  generation  of 
noise-induced  false  alarms  with  clutter-induced  false  alarms.  Arbitrary  values  were 
selected  in  order  to  highlight  the  behavior.  In  the  case  of  noise  (A)  the  false  alarms  are 
uniformily  distributed  in  range  while  for  clutter  (B)  the  occurrence  of  false  alarms 
rapidly  decreases  with  range. 


b 


FIGURE  8  -  An  example  of  generation  of  false  alarms  (Histogram) 
a)  Noise-induced,  b)  Sea-clutter-induced 
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8.0  CONCLUSION 

This  document  presented  a  model  for  the  simulation  of  a  surveillance  radar  for 
sensor  data  fusion  studies.  A  panoramic  pulse-MTI/Doppler  search  radar,  typical  of  the 
Canadian  Navy  surveillance  radars,  using  a  rotating  antenna  that  scans  360  degrees  in 
azimuth  was  considered.  The  proposed  model  is  more  flexible  and  more  generally 
applicable  than  the  models  published  so  far.  It  considers  a  wider  range  of  parameters 
which  affect  a  search  or  surveillance  radar’s  performance;  it  is  therefore  likely  to  yield 
more  “realistic”  results.  In  particular,  this  model  brings  the  following  improvements: 

•  an  extension  to  sea-clutter  models  to  account  for  ducting; 

•  better  identification  of  the  various  losses  for  a  typical  surveillance  radar; 

•  full  Parabolic  Equation  Method  (PEM)  computation  of  the  propagation  factor; 

•  a  new  model  to  generate  false  alarm  measurements  induced  by  sea  clutter; 

•  a  flexible  and  generic  model  to  represent  the  effects  of  Doppler  processing; 

•  a  simulation  architecture  allowing  refinement  in  modeling  as  well  as  real-time 

use  for  data  fusion  studies. 

Most  radar  simulations  have  exclusively  regarded  the  evaluation  of  target 
probability  of  detection  in  noise/clutter,  and  the  current  knowledge  of  the  clutter  process 
is  not  sufficiently  advanced  to  provide  accurate  representation  of  clutter-produced  false 
alarm  outcomes.  We  have  chosen  therefore  to  use  heuristic  methods  to  achieve  a  more 
likely  spreading  of  false  alarms  as  a  function  of  radar  characteristics  and  environment.  A 
probability  of  false  alarms,  P/ac  was  estimated  as  a  function  of  range  and  azimuth. 

A  two-level  processing  architecture  was  proposed  where  the  effects  of 
environment  and  Doppler  processing  are  precalculated  and  stored  in  tables:  rain  and  sea 
clutter  reflectivity  tables,  a  propagation  factor  table,  the  Doppler  processing  improvement 
factor  table  and  the  probability  of  false  alarms,  Pfac  table.  Tables  have  to  be  precomputed 
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for  the  environmental  conditions  and  radar  characteristics  for  the  cases  of  interest  to  be 
considered  by  CASE_ATTI. 

The  two-level  processing  architecture  provides  the  flexibility  of  using  various  non- 
standardized  pieces  of  software  to  precompute  sensitive  parameters  without  increasing 
complexity  and  therefore  affecting  the  execution  speed  of  the  CASE_ATTI.  Another 
advantage  is  that  these  external  modules  can  be  more  conveniently  upgraded  as  improved 
models  become  available.  The  model  was  particularly  suited  to  data  fusion  studies  where 
true  and  false  measurements  have  to  be  generated  while  the  scenario  is  running.  True 
measurements  are  generated  according  to  the  Pj  and  the  signal-to-interference  ratio  and 
false  measurements  are  generated  as  a  result  of  clutter  and  system  noise. 

This  radar  surveillance  model  is  an  important  step  towards  the  establishment  of  the 
CPF  sensor  suite  baseline  performance.  Two  radar  sensors  of  the  actual  CPF  sensor  suite 
can  be  represented  (SPS-49  and  Sea  Giraffe)  by  the  model  described  in  this  report.  This  is 
needed  to  support  the  ongoing  Multi-Sensor  Data  Fusion  (MSDF)  performance  evaluation 
study  in  the  CPF  context 
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